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Abstract 

Measurements on real traffic have revealed the existence of metastable states with 
very high flow. Such states have not been observed in the Nagel-Schreckenberg 
(NaSch) model which is the basic cellular automaton for the description of traffic. 
Here we propose a simple generalization of the NaSch model by introducing a velocity- 
dependent randomization. We investigate a special case which belongs to the so- 
called slow-to-start rules. It is shown that this model exhibits metastable states, 
thus sheding some light on the prerequisites for the occurance of hysteresis effects in 
the flow-density relation. 
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1 Introduction 



A few years ago, Nagel and Schreckenberg Q (NaSch) have proposed a probabilistic 
cellular automaton (CA) for the description of single-lane highway traffic^. Using 
very simple rules, this model is able to reproduce the basic phenomena encountered 
in real traffic, e.g. the occurance of phantom traffic jams ('jams out of nowhere'). 
The NaSch model is 'minimal' in the sense that every simplification of the rules no 
longer produces realistic results. On the other hand, for the description of more 
complex situations (e.g. multi-lane traffic, ramps)| or for a proper modelling of the 
'fine-structure' of traffic flow, additional rules have to be added and/or the basic 
rules have to be modified. 

The NaSch model |I| is a probabilistic cellular automaton. Space and time (and hence 
the velocities) are discrete. The road is divided into cells of length 7.5 m. Each cell 
can either be empty or occupied by just one car. The state of car j (J = 1, . . . , N) 
is characterised by its momentary velocity Vj (vj — 0,1, ... , v max ). The state of the 
system at time t+1 can be obtained from the state at time t by applying the following 
four rules to all cars at the same time (parallel dynamics): 

Rl: Acceleration: Vj(t) — ► Vj(t + ~) = min{u,(t) + l,v max } 

R2: Braking: if Vj(t + ~) > d,(t) then Vj(t + |) = dj(t) 

else Vj(t + Ti ) = Vj(t + I) 

R3: Randomization: Vj(t + §)—>■ Vj(t + 1) = max{0, Vj(t + |) — 1} 

with probability p 

R4: Driving: car j moves Vj(t + 1) cells. 

Here dj(t) denotes the number of empty cells in front of car j, i.e. the gap or headway. 
One timestep t — > t + 1 corresponds to approximately 1 sec in real time Jl|. 

In the spirit of modelling complex phenomena in statistical physics, the NaSch model 
does not try to describe traffic flow very accurately on a microscopic level. Macro- 
scopic effects observed in real traffic, e.g. the spontaneous formation of jams, can be 
understood by introducing just one simple stochastic parameter, the braking prob- 
ability p. Note that the motion of a single car might exhibit (on short timescales) 
unrealistic features, like stopping from maximum velocity within a few timesteps 
without any reason. However, after averaging over the motion of all cars or on long 
timescales, the NaSch model produces quite realistic results. Therefore one should 
not try to relate these large fluctuations to those observed in real traffic. 

Besides the CA models which are discrete in space and time, several other approaches 
to traffic flow have been discussed recently. Among these are space-continuous mod- 

1 For an overview of other approaches, see e.g. ||, [|, [j| 
2 For applications to urban traffic, see e.g. || fjj. 



2 



els in discrete time like the model of Krauss et. al. 0, §], as well as models continu- 
ous both in space and time, e.g. the macroscopic (fluid-dynamical) models |H]], the 



optimum- velocity model [11|, coupled- map models [|12[ and gas-kinetic models [13]. 
For further references we refer to [fj, [| |J. 

In the present paper we want to investigate hysteresis effects encountered in empiri- 
cal observations |L4|, [15], [16| . Such effects are related to the existence of metastable 
states in certain density regimes. In the NaSch model these states have not be ob- 
served. Here we want to present slightly modified models which are able to produce 
metastable states and hysteresis. As a consequence, the occurance of metastable 
states is therefore not related to the use of realistic braking rules, a continuum de- 
scription or deterministic models. 

In order to establish the existence of metastable states one can follow two basic 
strategies. In the first method, the density of cars is changed adiabatically by adding 
or removing vehicles from the stationary state at a certain density. Starting in the 
jamming phase (large densities) and removing cars, one obtains the lower branch 
of the hysteresis curve. On the other hand, by adding cars to a free flowing state 
(low densities), one obtains the upper branch. The second method does not require 
changing the density. Instead one starts from two different initial conditions, the 
megajam and the homogeneous state. The megajam consists of one large, compact 
cluster of standing cars. In the homogeneous state, cars are distributed equidistantly 
(with one large gap for incommensurate densities). 

In certain density regimes the fundamental diagram can consist of two branches. In 
the upper branch (with higher flow) there are almost no interactions between the 
cars and the system remains in a homogeneous, jam-free state. In the lower branch, 
however, the system is in a 'phase-separated' state, consisting of one large jam and 
a free-flowing part. 

Experimental observations |17| suggest that a reduction of the outflow from a jam 
compared to the maximum possible flowj] is an important ingredient for the occur- 
rance of metastable states with large lifetimes. Such a reduction can be implemented 
by so-called slow-to-start rules, where standing cars accelerate with lower probability 
than moving cars. It leads to a downstream flow with a relatively low density cor- 
responding to the lower branch of the hysteresis curve. Due to the reduction of the 
density in the outflow region of a jam, the flow in the jammed state is significantly 
lower than the flow of the homogeneous state at the same density. 

The paper is organized as follows: In Sect. H two models with slow-to-start rules 
are investigated. In Sect. |^ we introduce a slightly modified NaSch model with a 
velocity-dependent randomization. In this way, one incorporates the basic ingredients 
necessary to produce hysteresis. In the concluding Section ^ we summarize our results 
and compare with those for other models and real traffic. 



3 This reduction is often referred to as 'capacity drop'. 
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2 Models with slow-to-start rules 



In this Section we briefly present results for two models with slow-to-start (s2s) 
the T 2 model 



rules [IS 



19, 20 



and the BJH model pl| . These models have been 
introduced in order to model the restart behaviour of stopped cars in a more realistic 
fashion. As will be demonstrated below, these s2s rules are an important ingredient 
for the occurance of metastable states, although the authors of |19|, [2l|] did not realize 



the connection between s2s rules and metastability encountered in real traffic. 



2.1 T 2 model 

Takayasu and Takayasu (T 2 ) |19| have first suggested a CA model with a s2s rule. 



This rule has been generalized in |18[ as follows: A standing car with exactly one 
empty cell in front of it accelerates with probability q t = 1 — pt, while all other cars 
accelerate deterministically. The other update rules (R2-R4) of the NaSch model 
are unchanged, e.g. all cars are still subject to the randomization step. Due to this 
modification already for v max = 1 the particle-hole symmetry is broken. 

In Fig. |I] we show the fundamental diagram of the T 2 model with v max = 5, p = 0.01 
and pt = 0.75. The system size used for the simulation was L = 1000. In order to 
equilibrate the system, 10000 lattice updates have been performed. The data shown 
in Fig. p] represent an average over 100000 sweeps through the lattice. 




Figure 1: Fundamental diagram of the T 2 model (v n 



5, p = 0.01, p t = 0.75) 



obtained by starting from two different initial conditions, a completely jammed state 
(full line) and a homogeneous state (broken line). The peak of the full line is a 
finite-size effect. 



Comparing the simulation results with those for the NaSch model, two qualitative 
differences have to be discussed. First, the fundamental diagram has an inflection 
point in the high density regime and second the flow behaviour is non-unique in a 
density regime below the density of maximum flow p ma x- 
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The existence of an inflection point for large pt has first been noticed in [18| for the 
case v m ax = 1. In contrast, the fundamental diagram of the NaSch model is always 
convex. The different behaviour of the T 2 model for large pt and large densities p is 
related to the fact that space is not used as efficiently as in the NaSch model. One 
finds a similar behaviour in space-continuous models [^2j . There is some experimental 
evidence that in certain situations the shape of the fundamental diagram differs from 
the convex form. This behaviour of the average flow can be easily obtained tuning 



the parameter p t [0, [23|] . 

The non-unique behaviour of the flow for densities just below p max is due to the fact 
that for these densities the average flow still depends on the initial configuration. 
The measurements in Fig. [I] have been performed by applying the second method 
described in the introduction. The lower branch of the fundamental diagram corre- 
sponds to measurements starting from a initially completely jammed configuration 
while the upper branch has been obtained starting with a homogeneous initialization. 
It should be mentioned that the same result can also be obtained by applying the 
first method, i.e. by changing the density adiabatically. Since the flow depends on 
the history of the system, diagrams like those of Fig. [I] are usually called hysteresis 
curves. 

For low densities the stationary state consists of homogeneous configurations which 
are completely jam-free. At higher values of the global density the configurations 
contain one large jam. In contrast to the NaSch model, no spontaneous formation of 
jams in the outflow region of the large jam has been observed for the system sizes we 
took into account. Therefore the jammed states are phase separated states, unlike 
for the case of the NaSch model. 

It should be noted that the maximum value of the average flow still depends on 
the number of updates as well as on the system size. Therefore we can not ex- 
clude that the stationary value of the average flow is unique in the thermodynamic 
limit. Nevertherless the metastable states are extremely stable, even for large sys- 
tems. Moreover, in contrast to standard problems of statistical mechanics, not the 
thermodynamic limit is relevant for practical purposes, but the behavior of systems 
of finite length (note that 10000 lattice sites correspond to a road of length 75 km 
in reality). Since in reality roads and observation times are always finite, the above 
results are sufficient for all practical purposes. 

In Fig. we show the fundamental diagram for v max = 1, p = 0.5 and p t = 1, i.e. 
stopped cars can only move if there are at least two empty cells in front. Obviously 
completely blocked states exist for densities p > 0.5, where the number of empty cells 
in front of each car is smaller than two. Since fluctuations are absent in those states, 
they have an infinite lifetime. Therefore the flow in the stationary state is zero. In 
the region 0.5 < p < 0.66 states with a finite flow exist. Although these states are 
not stationary, one has to perform an extremely large number of update steps until 
the flow vanishes for large system sizes and densities slightly above p = 0.5, because 
the number of blocked configurations is very small compared to the total number of 
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configurations. Precisely at p = 0.5, the blocked state is unique and the typical time 
to reach this state diverges exponentially with the system size. Therefore we used 
very small systems in order to obtain the lower branch of the fundamental diagram. 




Figure 2: Fundamental diagram of the T 2 model (v Tl 



1, p — 0.5, p t = 1) for two 



system sizes. For densities slightly above p 
reached for the smaller system. 



0.5 the stationary state could only be 



Note that the mechanism for metastability in the case pt = 1 is different. The hys- 
teresis curve in Fig. ^ has been obtained by starting from two different homogeneous 
states, differing only in the velocity of the cars. The jammed branch is reached by 
starting with a configuration where all cars j have a velocity Vj(t = 0) = 0, whereas 
the upper branch corresponds to an initialization Vj(t = 0) = v max . In contrast, by 
following the first method described in the introduction, one would not find hystere- 
sis, since a megajam initialization would also yield the high-current branch. 

The reason for the occurance of hysteresis in the limit v max = 1 and pt — 1 is the 
existence of a "geometrical" phase transition. If the critical density is exceeded, the 
cars can no longer move since there is not enough free space. The "geometrical" 
phase transition makes it possible to find hysteresis even in the case v max = 1. 

Finally we want to remark that the p t = 1 limit of the T 2 model is some sense 
complementary to the cruise-control limit 24] of the NaSch model. In the T 2 model 
the completely blocked state is stabilized due to the absence of fluctuations whereas 
in the cruise-control limit one finds the absence of fluctuations for homogeneous states 
at low densities. 



2.2 BJH model 

The s2s rule of the T 2 model is a 'spatial' rule. The range of interaction for standing 
cars is larger than in the NaSch model and the restart behaviour depends only on the 
spatial arrangement of the vehicles. However, there are other ways of implementing 
a s2s behaviour. In the Benjamin- Johnson-Hui (BJH) model |2l| cars which had 
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to brake due to the next car ahead, will move on the next opportunity only with 
probability 1 — p s . Note that in contrast to the T 2 this slow-to-start rule requires 
'memory', i.e. it is a 'temporal' rule depending on the number of trials and not on 
the free space available in front of the car. 

For the BJH model no metastable states and hysteresis effects have been found until 
now since only the case v max = 1 has been investigated thoroughly [fjl, 



18 . For 



higher velocities we expect the occurance of metastable states, since also in the BJH 
model the outflow from a jam is smaller than the maximal flow. 




Figure 3: Fundamental diagram of the BJH model (v max = 5, p = 0.01, p s = 0.75) 
obtained using two different initial conditions, namely a completely jammed state 
(full line) and a homogeneous state (broken line). 



Our simulations show that for v max > 1 the overall behaviour of the BJH model is 
very similar to that of the T 2 model (for p t < 1). Therefore we are not going to 
discuss it here. In Fig. § we show a typical fundamental diagram. There is, however, 
no inflection point and the fundamental diagram is convex for all p s . This supports 
the view that the existence of an inflection point is a 'spatial' effect. 



3 NaSch model with velocity-dependent random- 
ization 

Here we present a simple generalization of the NaSch model which incorporates slow- 
to-start behaviour without introducing memory (like in the BJH model) or a longer- 
ranged interaction (like in the T 2 model). This new s2s rule is therefore neither 
temporal nor spatial. 

Instead, a velocity-dependent randomization (VDR) parameter p = p(v(t)) is intro- 
duced. This parameter has to be determined before the acceleration step Rl. For 
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simplicity we here study only the case 



p(v) 



Po v = 
p v > 



(i; 



which already contains the most important features of the general case |2q| . Since 



we are interested in hysteresis phenomena, we restrict ourselves to the case po < p. 
Note that for po = p the NaSch model is recovered. The cruise-control limit |24[ 
corresponds to the choice p(iWr) = and p{v ) = p for v < v max . 

In the following we will use a maximum velocity v max = 5, braking probability 
p = 1/64 of the moving cars and a higher value po = 0.75 for the braking probability 
of stopped cars. Simulation runs have been performed for periodic systems with 
L = 10000 lattice sites. 




Figure 4: Fundamental diagram of VDR model (v max — 5, pq — 0.75, p = 1/64, 
L = 10000). For comparison the fundamental diagrams of the NaSch model with 
p = 0.75 and p = 1/64 are given. 

Fig. |] shows the fundamental diagram of the modified model. Obviously the average 
flow J(p) can take two values in the density interval between p\ and p<i depending 
on the chosen initialization. The larger values of the average flow can be obtained 
using a homogeneous initialization of the system. The lower branch is obtained 
starting from a completely jammed state. Moreover, varying the particle number 
adiabatically, one can trace a hysteresis loop. One gets the upper branch by adding 
cars to the stationary state with p < pi and the lower one by removing cars from the 
stationary state with p > p 2 . For a fixed value of p, AJ = J{p2) — J(pi) depends 
linearly on po for wide a range of parameters. 

Increasing the system size, we observe a decrease of the density p 2 towards the branch- 
ing density p±. The jammed states become stable for densities p > p% even for global 
densities very close to p\. Again it should be noted that the homogeneous states 
have extremely long life-times and should therefore be relevant for realistic systems. 

It is instructive to compare the fundamental diagram of the VDR model with those 
of the corresponding NaSch models. For small densities p < 1 there are no slow 
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cars in the VDR model since interactions between cars are extremely rare. Here 
the flow is given by J(p) ~ p{v m ax — p), i-e. identical to the NaSch model with 
randomization p. For large densities 1 — p <C 1, on the other hand, the flow is given 
by J(p) ~ (1 — Po)(l ~ P) which corresponds to the NaSch model with randomization 
Po. For densities close to p = 1, only cars with velocities Vj — or Vj = 1 exist. The 
number of moving cars goes to zero so that asymptotically the flow is completely 
determined by po- 

The microscopic structure of the jammed states in the VDR model differs from 
those found in the NaSch model. While jammed states in the NaSch model contain 
clusters with an exponential size-distribution f23fl , one can find phase separation in 
the VDR model. The reason for this behaviour is the reduction of the outflow from 
a jam. If the outflow from a jam is maximal, any small jam in the free flow regime 
dissolves immediately since the outflow from such a jam is larger than the global 
flow. Therefore phase separation can not occur in that case. However, if the outflow 
from a jam is reduced, the density in the free flow regime is smaller than the density 
of maximum flow and cars can propagate freely in the low density part of the lattice. 
Due to the reduction of the density in the free flow regime, no spontaneous formation 
of jams is observable in the stationary state, if fluctuations in the free flow regime 
are rare. 

This picture is supported by a simple phenomenological approach. Obviously the 
flow in the homogeneous branch is given by Jhom = p(y max —p) = pvj, because every 
car can move with the free-flow velocity Vf. Assuming that the high density states 
are phase separated, we can obtain the second branch of the fundamental diagram. 
The phase separated states consist of a large jam and a free flow regime, where each 
car moves with velocity Vf. The density in the free flow regime pf is determined 
by the average waiting time T w = of the first car in the jam and Vf, because 
neglecting interactions between cars, the average distance of two consecutive cars is 
given by Ax = T w vj + 1 = pj . Using the normalisation L = Nj + N F Ax (Np(j) is 
the number of cars in the free flow regime (jam)) we find that the flow is given by 



Obviously pf is precisely the lower branching density pi, because for densities below 
Pf the jam-length is zero. It should be noted that this approach is only valid for 
Po 3> p and v max > 1. The condition p guarantees that interactions of cars due 
to velocity fluctuations are rare. As a consequence, the jam is compact in that limit. 
For increasing p, the jam becomes less dense. In the case v max = 1, cars can stop 
spontaneously, even in the free-flow regime. If po is large enough, these cars might 
initiate a jam. This is the basic reason why hysteresis is usually not observed for 



Measurements of the average flow show that the lower branch of the fundamental dia- 
gram is not stable near the density p±, if small system sizes are considered. Therefore 
we performed a more detailed stability analysis of the homogeneous and the jammed 




(2) 
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Figure 5: Time-dependence of the length Lj am (t) of the jam for p = 0.095. The left 
part of the figure shows the time evolution of the length L Jam (t) of one sample. The 
average (L Jam (t)) over 10000 samples (right part of the figure) shows an exponential 
decay. 



state near p\ and p 2 . Close to pi, the large jam present in the initial configuration 
dissolves and the average length (L Jam (t)) decays exponentially in time (Fig. |5|). It 
should be noted that this behaviour is not the consequence of a continuous "melting" 
of the large jam. In contrast, the jam-length Lj am (t) is strongly fluctuating without 
any obvious systematic time-dependence (Fig. ||). Once a homogeneous state without 
a jammed car is reached, no new jams are formed. Therefore the homogeneous state 
is stable near p 1 . For large system sizes the jammed states are stable for p > pi, 
also for densities only slightly above pi, because the average length of the jam is 
proportional to the system size, while the fluctuations grow sub-extensive. 




0.0 100.0 200.0 



Figure 6: Space-time diagram of the VDR model for p = 0.15, L = 400, p = 0.01 
and po = 0.5. The homogeneous initial state is not destroyed immediately, but after 
approximately 93000 lattice updates. In the outflow regime of the jam the density is 
reduced compared to the average density. 
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Analogous to the metastable jammed states near pi, homogeneous initializations for 
densities slightly above p 2 lead to metastable homogeneous states with short lifetimes. 
Fig. ^| shows the spontaneous formation of jams due to velocity fluctuations. The 
finite lifetimes of the homogeneous states are the qualitative difference between this 



model and the cruise-control limit p3 of the NaSch model, where the time evolution 



of homogeneous states at low densities is completely deterministic. 

For large system sizes the density difference Ap = p 2 — p\ decreases. This can be 
explained by looking at the mechanism of an emerging jam. Jams emerge due to 
velocity fluctuations in dense regions of the homogeneous states, where the distance 
dj + 1 between consecutive cars is less than v max . In these regions all following 
cars have to slow down if a car in front breaks in the randomization step. Over- 
reactions of following cars finally can cause jams. The probability to find clusters of 
an appropriate length for a given density is proportional to the system size. Therefore 
we expect that all homogeneous states are unstable for p > p\ in the thermodynamic 
limit. 

For higher values of p the lifetime of homogeneous states at densities p > p\ is 
very small. Nevertheless one can observe the same microscopic structure of the 
high density states as long as the outflow of a jam is sufficiently reduced. A rough 
estimate for the minimal difference between p and p necessary to observe phase 
separation can be obtained from the following arguments based on the features of 
a typical fundamental diagram of the VDR model at small p and large v max (see 
Fig. [Fp. In the following we assume that the fundamental diagram is perfectly linear 




Figure 7: Illustration of the flow branches a) in the presence of phase separation and 
b) without phase separation. 

up to p 2 , i.e. Jhom(p) = p(v max — p). As mentioned before, the phase-separated 
branch is described by J sep (p) = (1 — Po)(l — P) ( see Eq. (|])). In order to find phase 
separation, the flow in the phase-separated branch at density p 2 must be lower than 
the flow in the homogeneous branch, i.e. Jhamip^) > Jsepipz) (left part of Fig. |7|). 
For a situation as depicted in the right part of Fig. [7| (i.e. for Jhom{P2) < Jse P (p2)), 
no phase-separation would be found. In the following we approximate p 2 by p max , 



11 



the maximum of the corresponding NaSch model with randomization p. Using the 
estimate p max = (l—p)/(v max + l) f26j, the inequality for the current at p2 ~ p max can 
be used to obtain an estimate for the minimal difference po — P necessary to observe 
phase separation (see Fig. H). Near the deterministic limits (p = and p = 1) 
and for large values of v max , already a small difference p — p suffices to generate 
phase separation. For all nondeterministic cases the states of maximum flow are not 
completely homogeneous and the flow in the high density branch is somewhat lower 
than J sep (p) = (1 —po)(l — p) for larger values of p. In addition, the large jam is not 
compact anymore. 




Figure 8: Estimates for the minimum difference Po — p necessary to observe phase 
separation. 



In order to substantiate the picture developed in this section we also investigated 
correlation functions and the effect of perturbations (i.e. stopping a car temporarily 
to induce a jam) [p7|. The results are in full agreement with our interpretation and 



will be presented elsewhere [p8 



Another interesting result concerns the outflow from a jam. Our simulations show 
p8| that the outflow is (almost) independent of the density of cars and depends 



only on p. This in agreement with measurements on real traffic [|17]] where it was 
found that the outflow only depends on road and weather conditions and typical 
characteristics of the cars. Furthermore we have found that a structure consisting 
of two separated jams is rather stable and can exist for long times. This is also in 
agreement with experimental observations fllT . 

For v max = 1 we used the so-called car-oriented mean-field theory |29| in order to 
obtain an analytical description of the fundamental diagram. The results are in good 
agreement with simulations |p7 |. As expected, one finds no phase separation at any 
Po for p > 0. In the deterministic limit, however, a hysteresis curve can be observed 
for po > 0. The fundamental diagram constists of two linear branches and is in 
perfect agreement with simulations 
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4 Summary and Discussion 



The ability to describe the existence of hysteresis is a serious test for any traffic flow 
model. In this paper we have investigated several variants of the NaSch model and 
showed numerically the existence of metastable states in a density region close to the 
maximum flow. 

The hysteresis in the fundamental diagram is related to the existence of metastable 
states. The latter are a consequence of phase separation in a certain density regime. 
The reason for this phase separation, on the other hand, is a reduced outflow from 
jams which destabilizes clusters forming in the outflow region. Such a reduction can 
be incorporated most easily by introducing slow-to-start rules which try to model the 
restart behaviour of standing vehicles more realistically. The three models presented 
here used different s2s rules which can be classified as spatial, temporal and velocity- 
dependent, respectively. For maximum velocity v max > 1 the models exhibit a similar 
behaviour (see below). For v max = 1 the spatial s2s rule of the T2 model is exceptional 
since it leads to a existence of a phase transition to a completely jammed state in 
the limit p t — > 1 where p t is the s2s parameter. 

The NaSch model with velocity-dependent randomization shows the coexistence of 
phase separated and homogeneous states in a density interval near the density p 2 of 
maximum flow. Near p 2 interactions between cars become important and one can 
find spontaneous formation of jams. Therefore the reduction of the density in the 
outflow regime of a jam leads to stable phase separated states. The reduction of 
interactions between cars in the free flow regime can be confirmed by a phenomeno- 
logical approach, which gives very accurate results. In contrast to the cruise-control 



limit of the NaSch model EM , where also metastable states can be found, fluctuations 



are present in both coexisting states. 

The results presented here have interesting applications. In real traffic one is usually 
interested in stabilizing the homogeneous branch of the fundamental diagram in 
order to maximize the throughput. This can be done using signals to control the 
inflow, as in the case of the Lincoln tunnel in New York Results of simulations 



corresponding to such a situation will be presented in a future publication |2^ . 

In ]22| a family of space-continuous models has been investigated. Depending on 
the values of the maximum acceleration and deceleration three classes can be dis- 
tinguished. Class III shows no realistic behaviour since no spontaneous jams are 
formed. A continuous analogue of the NaSch model belongs to class II (high decla- 
ration limit). Here one finds spontaneous formation of jams, but no metastable states. 
Finally, in class I spontaneous jamming as well as metastable states are found. On 
a macroscopic level, models in classes I and II can be distinguished by the ordering 
of the densities p/, the density of the outflow from a jam, and p c , where the density 
of a homogeneous flow becomes unstable^. The discrete models discussed here have 

4 Oftcn it is not easy to define these quantities properly, especially in the presence of noise. 
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very simple braking rules, but nevertheless exhibit the behaviour of class III. The 
existence of metastable states is therefore not related to the use of continuous space 
coordinates or "realistic" braking rules || [J (e.g. anticipation of the behaviour of 
the driver ahead). 

Already the NaSch model has a tendency towards the formation of metastable states. 
The outflow from a megajam is maximal only in the deterministic case p = and 
in the so-called cruise-control limit |[24|| . In the generic case p > the outflow is 



not maximal [31|, but organizes itself towards the 'critical' current J c where the 



c 



correlation length becomes maximal f2"B| . Since the difference between J m ax and J, 
is rather small, it is very difficult to observe metastable states in the NaSch model. 

The considerations in this paper show the flexibility of the CA approach to traffic 
flow problems. A rather simple and natural extension of the rules of the NaSch 
model allows us to describe the formation of metastable states in the fundamental 
diagram. The introduction of a velocity-dependent randomization p(v) makes it 
possible to control the properties of the free flow and congested flow independently. 
Experimentally it has been found [17] that for real traffic the reduction of the outflow 
J ou t compared to the maximum flow J ma x is approximately J ma x/Jout ~ 1-5. This 
value can be used to determine a realistic value of p — p. 

If one is willing to give up an overly realistic description of the interactions between 
the vehicles one can obtain rather simple CA models capable of describing even the 
fine-structure of traffic flow in a satisfactory way. 
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